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Abstract 



In systems with a mixed phase space, where regular and chaotic motion coexists, 
regular states are coupled to the chaotic region by dynamical tunneling. We give an 
overview on the determination of direct regular-to-chaotic tunneling rates using the 
fictitious integrable system approach. This approach is applied to different kicked 
Q ' systems, including the standard map, and successfully compared with numerical 

^ ■ data. 

£h ! This text corresponds to Chapter 6 of the book: Dynamical Tunneling - Theory and 

Experiment, edited by S. Keshavamurthy and P. Schlagheck [Taylor and Francis 
CRC (2011)] J7]/. For a more extensive exposition see [Phys. Rev. E 82, 056208 
(2010); \arXiv:1009.0\~T8\ 2]. 
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Tunneling is one of the most fundamental manifestations of quantum mechanics. For ID 
O . systems the theoretical description is well established by the Wentzel-Kramers-Brillouin 
(WKB) method and related approaches [21 E]. However, for higher dimensional systems 
no such simple description exists. In these systems typically regular and chaotic motion 
coexists and in the two-dimensional case regular tori are absolute barriers to the motion. 
Quantum mechanically, the eigenfunctions either concentrate within the regular islands 
or in the chaotic sea, as expected from the semiclassical eigenfunction hypothesis [H El 
|6]. These eigenfunctions are coupled by so-called dynamical tunneling [7] through the 
dynamically generated barriers in phase space. In particular, this leads to a substantial 
enhancement of tunneling rates between phase-space regions of regular motion due to the 
presence of chaotic motion, which was termed chaos-assisted tunneling [HI EJ EE] • Such 
dynamical tunneling processes are ubiquitous in molecular physics and were realized with 
microwave cavities [UJ and cold atoms in periodically modulated optical lattices [T2| IT3]. 

In the quantum regime, h e Q < A, where the effective Planck constant h e R is smaller 
but still comparable to the area of the regular island A, the process of tunneling into the 
chaotic region is dominated by a direct regular-to- chaotic tunneling mechanism [IU E3 EEl 
ITT] . For the prediction of tunneling rates the fictitious integrable system approach was in- 
troduced recently [17] . It relies on a fictitious integrable system [HI [15] that resembles the 
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regular dynamics within the island under consideration. The approach has been applied 
to quantum maps [17], billiard systems [IB], and the annular microcavity [TH]. In the semi- 
classical regime, h e g <C A, however, the direct tunneling contribution alone is typically 
not sufficient to describe the observed tunneling rates, because nonlinear resonance chains 
within the regular island lead to resonance-assisted tunneling [201 EH E21 [23] . Recently, a 
combined prediction of dynamical tunneling rates from regular to chaotic phase-space re- 
gions was derived [21] , which combines the direct regular-to-chaotic tunneling mechanism 
in the quantum regime with an improved resonance- assisted tunneling theory in the semi- 
classical regime. In this text we concentrate on the fictitious integrable system approach 
for the theoretical description of the direct regular-to-chaotic tunneling mechanism and 
how it can be applied analytically, semiclassically, and numerically to quantum maps. 

This text is organized as follows: We start by defining different convenient kicked 
systems, their quantization, and the numerical determination of tunneling rates. In Sec. [3] 
we describe the approach to determine dynamical tunneling rates by using a fictitious 
integrable system. In Sec. H] we apply this theory to the previously introduced systems 
and compare its results with numerical data. We conclude with a brief summary. 



Kicked systems are particularly suited to study classical and quantum effects appearing in 
a mixed phase space as they can be easily treated analytically and numerically. In contrast 
to time-independent systems, where at least two degrees of freedom are necessary to break 
integrability, one-dimensional driven systems can show chaotic motion. We consider time- 
periodic one- dimensional kicked systems 



which are described by the kinetic energy T(p) and the potential V(q) which is applied 
once per kick period r = 1. The classical dynamics of a kicked system is given by its 
stroboscopic mapping, e.g., evaluated just after each kick 



It maps the phase-space coordinates after the nth kick to those after the (n + l)th kick. 
For the stroboscopic mapping we consider a compact phase space with periodic boundary 
conditions for q G [—1/2, 1/2] and p G [—1/2, 1/2]. The phase space generally consists 
of regions with regular motion surrounded by chaotic dynamics, see Fig. D^b). We focus 
on the situation of just one regular island embedded in the chaotic sea, Fig. |2]^a),(b). At 
the center of the island one has an elliptic fixed point, which is surrounded by invariant 
regular tori. Classically, a particle within a regular island will never enter into the chaotic 
region and vice versa, i.e. tori form absolute barriers to the motion. 



2 Kicked systems 




(i) 



q n+ i = q n + T'ipn), 

Pn+l =Pn~ V'(q n+1 ). 



(2) 
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Figure 1: Phase space of the standard map for (a) k = 0.3, (b) k = 2.4, and (c) k = 8. 
In (a) the system is nearly integrable. Most of the phase space is filled by invariant tori 
(lines). The second figure (b) shows a mixed phase space with regular motion (lines) 
and chaotic motion (dots). The four small islands surrounding the central island form a 
nonlinear resonance chain. In the chaotic sea partial barriers exist, which is reflected by 
the different densities of chaotic points. For large values of k as in (c) the phase space 
seems macroscopically chaotic. 



2.1 Standard map 

The paradigmatic model of an area preserving map is the standard map, defined by Eq. (j2J) 
with the functions 

T'{p) = p, (3) 
V'(q) = |^sin(27rg). (4) 

For k = the system is integrable and the dynamics takes place on invariant tori with 
constant momentum. For small k > 0, see Fig. QJa), the system is no longer integrable. 
However, as stated by the Kolmogorov- Arnold- Moser (KAM) theorem [25] , many invariant 
tori persist. In between these tori sequences of stable and unstable fixed points emerge, 
which is described by the Poincare-Birkhoff theorem [25]. The stable fixed points of such 
a sequence form so-called nonlinear resonance chains, see Fig. H(b), while in the vicinity 
of the unstable fixed points typically chaotic dynamics is observed. For larger k more 
and more regular tori break up and the chaotic dynamics occupies a larger area of phase 
space. At k ~ 3 one large regular region ("regular island") is surrounded by a region 
of chaotic dynamics ("chaotic sea"). At the border of the regular island a hierarchical 
region exists, where self-similar structures of regular and chaotic dynamics are found on 
all scales. Also in the chaotic sea additional structures can be found. Here, partial barriers 
limit the classical flux between different regions of chaotic motion caused by the stable 
and unstable manifolds of unstable periodic orbits or caused by cantori, which are the 
remnants of regular tori [27]. For even larger k the regular islands disappear and the 
motion seems macroscopically chaotic, see Fig. QJc). 

All these structures in phase space, such as nonlinear resonances, the hierarchical 
regular-to-chaotic transition region, and partial barriers within the chaotic sea, have an 
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influence on the dynamical tunneling process from the regular island to the chaotic sea. 
In order to quantitatively understand dynamical tunneling in systems with a mixed phase 
space, it is helpful to first consider simpler model systems. By designing maps such that 
the phase-space structures can be switched on and off one-by-one it is possible to study 
their influence on the dynamical tunneling process separately. Afterwards the results can 
be tested using the generic example of the standard map. 



2.2 Designed maps 

Our aim is to introduce kicked systems which are designed such that their phase space is 
particularly simple. It shows one regular island embedded in the chaotic sea, with very 
small nonlinear resonance chains within the regular island, a narrow hierarchical region, 
and without relevant partial barriers in the chaotic component. For such a system it is 
possible to study the direct regular-to-chaotic tunneling process without additional effects 
caused by these structures as long as h e s is big enough. To this end we define the family 
of maps T>, according to Eq. (j2J), with an appropriate choice of the functions T'(p) and 
V'(q) [2E1 EH EH 123. For this we first introduce 

t'(p) = ^±(l-2p) for 0<±p<^ (5) 
v'(q) = —rq + Rq 2 for — - < q < — (6) 

with < r < 2 and R > 0. Considering periodic boundary conditions the functions t'(p) 
and v'(q) show discontinuities at p — 0, ±1/2 and q = ±1/2, respectively. In order to 
avoid these discontinuities we smooth the periodically extended functions v'(q) and t'(q) 
with a Gaussian, 

resulting in analytic functions 

T'{p) = J dzf(z)G(p-z), (8) 

V\q) = J dzv'(z)G(q-z), (9) 

which are periodic with respect to the phase-space unit cell. With this we obtain the maps 
T> depending on the parameters r, R, and the smoothing strength e. The smoothing e 
determines the size of the hierarchical region at the border of the regular island. Tuning 
the parameters r and R one can find situations, where all nonlinear resonance chains 
inside the regular island are small. 

For R = both functions v'(q) and t'(p) are linear in q and p, respectively. In this case 
we find a harmonic oscillator-like regular island with elliptic invariant tori and constant 
rotation number. We choose the parameters r = 0.46, R = 0, and e = 0.005 for which 
the phase space of the resulting map is shown in Fig. E£a)- 

In typical systems the rotation number of regular tori changes from the center of the 
regular region to its border which typically has a non-elliptic shape. Such a situation can 
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Figure 2: Phase space of (a) the designed map 2\ with a harmonic oscillator-like regular 
island, (b) V d with a deformed regular island, and (c) V rs with a regular stripe. For 
parameters see text. 



be achieved for the family of maps T> with the parameter R ^ 0. For most combinations 
of the parameters r and R resonance structures appear inside the regular island. They 
limit the /i e ff-regime in which the direct regular-to-chaotic tunneling process dominates. 
Hence, we choose a situation in which the nonlinear resonances are small such that their 
influence on the tunneling process is expected only at small h e s. The phase space of the 
map V d , obtained for r = 0.26, R = 0.4, and e = 0.005, is shown in Fig. [2](b). 

Another designed kicked system was introduced in Refs. [5UI EU E2] • Here the regular 
region consists of a stripe in phase space, see Fig. Ulc). In our notation the mapping, 
Eq. ([2]), is specified by the functions 



V'(q) 



T'ip) 



— ( 8iraq + di — d 2 + -[87rag — u + di] tanh[6(87rg — qa)] 
2tt \ 2 



-8iiaq + oj + d 2 ] tanh[6(87rg + q d )} 



H 



■— sin(27rp). 

07T 



(10) 

fill 



The kinetic energy T(p) is periodic with respect to the phase space unit cell. We label the 
resulting map with a regular stripe by V rs using the parameters a = 5, b = 100, d\ = —24, 
d 2 = —26, u = 1, qd = 5, and n = 3. 

The map T> TS is similar to the system Dho as it also destroys the integrable region 
by smoothly changing the function V'(q), here at \q\ = qd/(Sn). For \q\ < qd/{8n) the 
potential term is almost linear while it tends to the standard map for \q\ > qd/(8n). The 
parameter b determines the width of the transition region. 
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2.3 Quantization 



Quantum mechanically kicked systems are described by the unitary time-evolution oper- 
ator U t331. It can be written as U = UyUx with 



-iV(q)/h eB 



-iT(p)/h B e 



U v = e 
U T = e 

Its eigenstates \ip n ) and quasi-energies ip n are determined by 



U\lf> n ) 



(12) 
(13) 

(14) 



We consider a compact phase space with periodic boundary conditions in q and p. This 
implies that the effective Planck constant can take only the discrete values 



(15) 



where N is the dimension of the Hilbert space. In addition the position and momentum 
coordinates are quantized according to 



Qk 
Pi 



k 1 
N~2 : 

L _ I 

N ~~ 2 



(16) 
(17) 



with k,l = 0, . . . ,JV- 1. 

For systems with a mixed phase space, in the semiclassical limit (/i cff — > 0), the 
semiclassical eigenfunction hypothesis jH El |6] implies that the eigenstates \ip n ) can be 
classified as either regular or chaotic according to the phase-space region on which they 
concentrate, see Fig. |3j In order to understand the behavior of eigenstates away from 
the semiclassical limit, i.e. at finite values of the effective Planck constant h c g, one has to 
compare the size of phase-space structures with h e g. 

The so-called regular states are concentrated on tori within the regular island and 
fulfill the Bohr- Sommerf eld type quantization condition 



p dq = h eS [m 
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Figure 3: Husimi representation of eigenfunctions of the map for h c g = 1/50. The 
regular ground state m = and the first three excited states m = 1,2,3 as well as a 
chaotic state are shown. 
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For a given value of h eS there exist N teg of such regular states, where N TCg = [A/ h cS + 1/2 J 
and A is the area of the regular island. The = N — N Teg chaotic states extend over the 
chaotic sea. Due to dynamical tunneling, however, the regular and chaotic eigenfunctions 
of U always have a small component in the other region of phase space, respectively. This 
is most clearly seen for hybrid states which even have the same weight in each of the 
components as they are involved in an avoided level crossing. 

2.4 Numerical determination of tunneling rates 

The structure of the considered phase space, with one regular island surrounded by the 
chaotic sea, allows for the determination of tunneling rates by introducing absorption 
somewhere in the chaotic region of phase space. For quantum maps, this can for example 
be realized by using a non-unitary open quantum map [35] 

U° = PUP, (19) 

where P is a projection operator onto the complement of the absorbing region. 

While the eigenvalues of U are located on the unit circle the eigenvalues of U° are 
inside the unit circle as U° is sub-unitary. The eigenequation of U° reads 

u°m = Zn m (20) 

with eigenvalues 




The decay rates are characterized by the imaginary part of the quasi-energies in Eq. f!2ip 
and one has 

lm = -2 log \z m \ « 2(1-1^1). (22) 

If the chaotic region does not contain partial barriers and shows no dynamical local- 
ization, it is justified to assume that the rate of escaping the regular island is equal to 
the rate of leaving through the absorbing regions located in the chaotic sea. Then, the 
location of the absorbing regions in the chaotic part of phase space has no effect on the 
tunneling rates. 

In generic systems, however, partial barriers will appear in the chaotic region of phase 
space. The additional transition through these structures further limits the quantum 
transport such that the calculated decay through the absorbing region occurs slower than 
the decay from the regular island to the neighboring chaotic sea. Similarly dynamical 
localization in the chaotic region may slow down the decay. The influence of partial 
barriers and dynamical localization on the tunneling process is an open problem. Here 
we will suppress their influence, if necessary, by moving the absorbing regions closer to 
the regular island. 

In Fig. H] tunneling rates are determined numerically by opening the system for the 
map 2\o (dots). These agree very well with an analytical prediction (lines), Eq. (1271) . 
which is derived in the next section. 
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Figure 4: Numerical data (dots) for m < 8 for the map with a harmonic oscillator-like 
island. Comparison with Eq. fl27|) (lines), see Sec. 13.11 for its derivation. The insets show 
Husimi functions of the regular states for the quantum numbers m = and m = 5 at 
l/h eS = 30. 

3 Tunneling using a fictitious integrable system 

Dynamical tunneling in systems with a mixed phase space couples the regular island 
and the chaotic sea, which are classically separated. This coupling can be quantified by 
tunneling rates 7 m which describe the decay of regular states to the chaotic sea. To define 
these tunneling rates one can consider a wave packet started on the mth quantized torus 
in the regular island coupled to a continuum of chaotic states, as in the case for an infinite 
chaotic sea or in the presence of an absorbing region somewhere in the chaotic sea. Its 
decay e _7m * is described by a tunneling rate j m . For systems with a finite phase space 
this exponential decay occurs at most up to the Heisenberg time th = h eS /A ch , where 
A c h is the mean level spacing of the chaotic states. 

In the quantum regime, h e g < A, where h e g is smaller but comparable to the area A of 
the regular island, the rates 7 m are dominated by the direct regular-to-chaotic tunneling 
mechanism, while contributions from resonance-assisted tunneling are negligible. We 
concentrate on situations where additional phase-space structures within the chaotic sea 
are not relevant for tunneling. In the following we derive a prediction for the direct 
regular-to-chaotic tunneling rates using the fictitious integrable system approach [Tf] , 
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3.1 Derivation 



In order to find a prediction for the direct regular-to-chaotic tunneling rates we decompose 
the Hilbert space of the quantum map U into two parts which correspond to the regular 
and chaotic regions. While classically such a decomposition is unique (neglecting tiny 
phase-space structures), quantum mechanically this is not the case due to the uncertainty 
principle. We find a decomposition by introducing a fictitious integrable system U Teg 
(a related idea was presented in Refs. [HI [E2 EE])- It has to be chosen such that its 
dynamics resembles the classical motion corresponding to U within the regular island 
as closely as possible and continues this regular dynamics beyond the regular island of 
U, see Fig. [5j(a),(d). The eigenstates |^reg) °f ^reg are purely regular in the sense that 
they are localized on the mth quantized torus of the regular region and continue to decay 
beyond this regular region, see Fig. [3(e). This is the decisive property of |^™ g ) which 
have no chaotic admixture, in contrast to the predominantly regular eigenstates of U, see 
Fig. [5j(b). The explicit construction of U reg is discussed in Sec. 13.21 

With the eigenstates \i/)™ g ) of U reg , U reg \ip™ g ) = e'^l^^g), we define a projection 
operator 

7V reg -l 

Pr Cg := J2 IC g ><C g l, (23) 

m=0 

using the first A^ rcg regular states of U rcg which approximately projects onto the regular 
island corresponding to U. The orthogonal projector 

P ch := 1 - P rcg (24) 

approximately projects onto the chaotic phase-space region. These projectors, P reg and 
P c h, define our decomposition of the Hilbert space into a regular and a chaotic subspace. 

Introducing an orthonormal basis \ip c h) i n the chaotic subspace we can write P c h = 
S c h lV , ch)('0ch|- Here we sum over all N ch = N — N rcg states (V'ch), see Fig. [5j(f) for an 
example. Hence, the coupling matrix element v c ^ m between a purely regular state \ip™ g ) 
and any chaotic basis state \ijjch) is 

V ch , m = (^ ch |£/|Cg>- (25) 

The tunneling rate is obtained using a dimensionless version of Fermi's golden rule, 

7m = Kh, m | 2 , (26) 

ch 

where the sum is over all chaotic basis states l^ch) and thus averages the modulus squared 
of the fluctuating matrix elements v c ^ m . Here we apply Fermi's golden rule in the case of a 
discrete spectrum, which is possible if one considers the decay e~ lmt up to the Heisenberg 
time th = /i c ff/A c h only. Inserting Eq. (l2"5j) we obtain 

l m = ||Pht/|Cg)H 2 = 11(1 - Pre S W™ g )\\ 2 (27) 

as the basis of all our following investigations. It allows for the prediction of tunneling rates 
from a regular state localized on the mth quantized torus to the chaotic sea. Equation (1271) 
confirms the intuition that the tunneling rates are determined by the amount of probability 



9 



u 





reg/ 




Figure 5: (a) Illustration of the classical phase space corresponding to some mixed quan- 
tum map U together with the Husimi representation of (b) a regular and (c) a chaotic 
eigenstate of U which both have a small component in the other region, (d) Illustration 
of the classical phase space of the fictitious integrable system U Teg . (e) Eigenstates |"0 re g) 
of U reg are purely regular while (f ) the basis states are localized in the chaotic region 
of phase space. 



of IV'reg) ^ na ^ * s transferred to the chaotic region after one application of the time evolution 
operator U . In Eq. (12?]) properties of the fictitious integrable system U reg and the chaotic 
projector P c h = 1 — P reg enter, which rely on the chosen decomposition of Hilbert space. 

In cases where one finds a fictitious integrable system U ieg which resembles the dynam- 
ics within the regular island of U with very high accuracy, Eq. (T27J can be approximated 
as 

1m « ||(tf-^eg)IC g )ir, (28) 

using P reg U\ip™ g ) ps P TCg U TCg \ip™ g ) = U Teg \ip™ g ) . Instead of the operator product P ch U in 
Eq. ( 127)1 the difference U — U rcg enters in Eq. (1251) . It allows for further derivations, which 
are presented in Sec. 13.31 

3.2 Fictitious integrable system 

The most difficult step in the application of Eq. ( |27j) to a given system is the determination 
of the fictitious integrable system U reg . On the one hand its dynamics should resemble 
the classical motion of the considered mixed system within the regular island as closely 
as possible. As a result the contour lines of the corresponding integrable Hamiltonian 
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ifreg, Fig. Efd), approximate the KAM-curves of the classically mixed system, Fig. Efa), 
in phase space. This resemblance is not possible with arbitrary precision as the integrable 
approximation does not contain, e.g., nonlinear resonance chains and small embedded 
chaotic regions. Moreover, it cannot account for the hierarchical regular-to-chaotic tran- 
sition region at the border of the regular island. Similar problems appear for the analytic 
continuation of a regular torus into complex space due to the existence of natural bound- 
aries [3S1 E7J EQ1 ED EH I2D1 12D 122] • However, for not too small h cS , where these small 
structures are not yet resolved, an integrable approximation with finite accuracy turns 
out to be sufficient for a prediction of the tunneling rates. 

On the other hand the integrable dynamics of H reg should extrapolate smoothly be- 
yond the regular island of H. This is essential for the quantum eigenstates of H Icg to 
have correctly decaying tunneling tails which are according to Eq. (1271) relevant for the 
determination of the tunneling rates. While typically tunneling from the regular island 
occurs to regions within the chaotic sea close to the border of the regular island, there 
exist other cases, where it occurs to regions deeper inside the chaotic studied 
in Ref. [16]. Here H reg has to be constructed such that its eigenstates have the correct 
tunneling tails up to this region. 

For quantum maps we determine the fictitious integrable system in the following way: 
We employ classical methods, see below, to obtain a one-dimensional time-independent 
Hamiltonian H reg (q,p) which is integrable by definition and resembles the classically reg- 
ular motion corresponding to the mixed system. After its quantization we obtain the 
regular quantum map U reg = e~ lHrcs ^ hcS which has the same eigenf unctions l^g) as -^reg- 
For the numerical evaluation of Eq. (127)) we use P ch = 1 — P TCg = 1 — l^regX^regl; 
according to Eqs. (1231) and (l2~3| . where the sum extends over m = 0, 1, . . . , iV re g — 1. 

Two examples for the explicit construction of H rcg will be mentioned below. Note, that 
also other methods, e.g. based on the normal-form analysis [391 HQ] or on the Campbell- 
Baker-Hausdorff formula [JT] can be employed in order to find H rcg . For the example 
systems considered in this manuscript, however, they show less good agreement. 

One possible choice for the determination of the fictitious integrable system for quan- 
tum maps is the Lie-transformation method [12]. It determines a classical Hamilton 
function as a series, 

K 

i=i 

see Ref. [21] and Fig. [6]^a) for an example. Typically, the order of the expansion K can 
be increased up to 20 within reasonable numerical effort. The Lie-transformation method 
provides a regular approximation H Teg which interpolates the dynamics inside the regular 
region and gives a smooth continuation into the chaotic sea. At some order K the series 
should diverge due to the nonlinear resonances inside the regular island. For strongly 
perturbed systems, such as the standard map at n > 2.5, the Lie-transformation method 
may not be able to reproduce the regular dynamics of U, see Fig. Mjo). Here we use a 
method [TTJ based on the frequency map analysis [13]. 

An important question is whether the direct tunneling rates obtained using Eq. ()2"7j) 
depend on the actual choice of H rcg and how these results converge depending on the 
order K of its perturbation series. Ideally one would like to use classical measures, which 
describe the deviations of the regular system H mg from the originally mixed one, to 
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Figure 6: Application of the Lie transformation method: (a) Orbits (gray) of the map 
an d of the corresponding integrable system (thin lines) of order K = 15. Here H Teg 
accurately resembles the regular dynamics of U. (b) Orbits of the standard map for 
k = 2.9 (gray) and of the corresponding integrable system (thin lines) of order K = 7. 
Here H veg does not accurately resemble the regular dynamics of H. 



predict the error of Eq. (I27j) for the tunneling rates. However, these classical measures 
can only account for the deviations within the regular region but not for the quality of 
the continuation of H reg beyond the regular island of U. Currently, the quality of an 
integrable system can be estimated a posteriori by comparison of the predicted tunneling 
rates with numerical data. It remains an open question how to obtain a direct connection 
between the error on the classical side and the one for the tunneling rates. 

3.3 Approximate fictitious integrable system 

For an analytical evaluation of the result, Eq. ( 12 7p . of the direct regular-to-chaotic tun- 
neling rates we approximate the fictitious integrable system U rcg by a kicked system, 
U TCg = UyUr or U TCg = UyUf with 

U v = e-^/*"*, (30) 
U f = e - if(p)/ ^«. (31) 

Here the functions V(q) and T(p) are a low order Taylor expansion of V(q) and T(p), 
respectively, around the center of the regular island. Note, that the classical dynamics 
corresponding to U TCg is typically not completely regular. Still the following derivation 
is applicable if U Teg has the following properties: (i) Within the regular island it has 
an almost identical classical dynamics as U, including nonlinear resonances and small 
embedded chaotic regions, (ii) It shows mainly regular dynamics for a sufficiently wide 
region beyond the border of the regular island of U. 

Now we consider the specific case U reg = U v Ut and assume that both properties (i) 
and (ii) are fulfilled. As the dynamics of U and U Teg are almost identical within the regular 
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island of U, the approximate result, Eq. (|28|) . can be applied with U reg replaced by U XI 
giving 



1m « ||(C/-^T)|^ 

= n(c/44-i)c/^ T ic g )ir- 



(32) 
(33) 



We now use that l^^g)' which is an eigenstate of the exact U reg and H veg , is an approximate 
eigenstate of U v Ut, leading to UyUxl^g) ~ e 1¥, "g|^ g ). With this we obtain 



/m. ~ 

In position representation this reads 

N-l 



2 x; 



fc=0 



1 — cos 



off 



(34) 



(35) 



where AV^(g) := V^(g) — V(g) and = k/N — 1/2. In the semiclassical limit the sum in 
Eq. (|3"5|) can be replaced by an integral 



1/2 



7n 



dg|C g (<?) 



1 — COS 



-1/2 



/AV%)\ 

V ^cff / 



(36) 



where we integrate over the whole position space q G [—1/2, 1/2]. Note, that for the 
complementary situation, where U rcg ~ UyUf is used in Eq. (j2"81) . a similar result can be 
obtained in momentum representation 



1/2 



4pK(p)| ! 



1 — cos 



-1/2 



fAT(p)Y 



(37) 



with AT(p) := T(p) - T(p). 

We now use a WKB expression for the regular states (VCg)- F° r simplicity we restrict 
to the case H Teg = p 2 /2 + W(q) leading to 



<c(ff) 



Co' 



27r|p(g)| 



exp 



1 



b(g')|dg' 



(3« 



which is valid for q > q r m . Here q r m is the right classical turning point of the mth quantizing 
torus, u is the oscillation frequency, and p(q) = ^j2(E™ g — W(q)). The eigenstates ip™ g (q) 
decay exponentially beyond the classical turning point q r m . The difference of the potential 
energies AV(q) = V(q) — V(q) approximately vanishes within the regular region and 
increases beyond its border to the chaotic sea. Hence, the most important contribution in 
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Eq. ( )36|) arises near the left or the right border, q\ or ql, of the regular island. For q > ql 
we rewrite the regular states 



/ 



lC(?) « C K (% r )exp 



C g (??) exp 



\ 



fteff J 



<ii; 



J 



I pig) 

p(q) 



h, 



eft 



-(<?-<?Db(<?DI 



(39) 
(40) 



where in the last step we use p(q) ~ p(ql) in the vicinity of the border. 

In order to use Eq. fl36|) we split the integration interval in two parts, such that 
lm = 1m + 7mi corresponding to the contributions from the left and the right. We now 
approximate 



AV(q) 



with some constant c b and find 



, q r m <q<ql 

Cb(q ~ ql) , q> q r b 



(41) 



2MC g (% r )l 2 



[1 — cos(q>x)] dx 



7T 



mM)[ 



(42) 



(43) 



where x = (q — ql)/h e g and x max = (1/2 — ql)/h e g. In the semiclassical limit x max — > oo 
and for fixed quantum number m the integral in Eq. (|42|) . 



-2^IK?PI[i _cos(c 6 x)]dx, 



(44) 



becomes an h e ^- independent constant. The tunneling rate 7^ is proportional to the square 
of the modulus of the regular wave function at the right border ql of the regular island. 
With Eq. (JHHD we obtain 
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eft 



27r 2 |p(g j 



exp 



V 



2 



\ 

\p{q')\dq' 

) 



(45) 



A similar equation holds for 7^. 

As an example for the explicit evaluation of Eq. (1451) we now consider the harmonic 
oscillator H rcg (p, q) = p 2 /2 + u 2 q 2 /2, where u denotes the oscillation frequency and gives 
the ratio of the two half axes of the elliptic invariant tori. Its classical turning points 
q^j = ±y/2E m /u, the eigenenergies E m = h c guj{m + 1/2), and the momentum p(q) = 
a/ 2E m — q 2 oo 2 are explicitly given. Using these expressions in Eq. ( 14"5]) and j m = 2Y m we 
obtain 



In 




(3 m - a m In 



(46) 
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as the semiclassical prediction for the tunneling rate of the mth regular state, where 
a m = (m + l/2)(A/h cS )~ 1 , j3 m = y/1 — a m , and A is the area of the regular island. The 
exponent in Eq. (146 p was also derived in Ref. [H]. The prefactor 

(47) 

can be estimated semiclassically by solving the integral, Eq. OH]), for x max — > oo. For a 
fixed classical torus of energy E one obtains 




2\p(q r b )\ + 4 



2 ■ 



(48) 

With this prefactor the prediction Eq. f H6|) gives excellent agreement with numerically 
determined data over 10 orders of magnitude in 7, see Fig. [7J For a fixed quantum number 
m the energy E m goes to zero in the semiclassical limit such that one can approximate 
|p(q£)| ~ uql in Eq. f]4"8"]) which does not depend on h cS . 

Let us make the following remarks concerning Eq. (j4"6"]) : The only information about 
this non-generic island with constant rotation number is A/h e g as in Ref. [15]. In contrast 
to Eq. ( l2Tj) it does not require further quantum information such as the quantum map U. 
While the term in square brackets semiclassically approaches one, it is relevant for large 
h c g. In contrast to Eq. (13T)]) . where the chaotic properties are contained in the difference 
V(q) — V(q), they now appear in the prefactor c via the linear approximation of this 
difference. 

In the semiclassical limit the tunneling rates predicted by Eq. ( 146]) decrease exponen- 
tially. For h c ff — > the values a m go to zero and f3 m to one, such that 7 ~ Q~ 2A / h ^ rema i n s 
which reproduces the qualitative prediction obtained in Ref. [14] . We find that the non- 
universal constant in the exponent is 2 which is comparable to the prefactor 3 — In 4 f=s 1.61 
derived in Refs. [T5] US] . However, our result shows more accurate agreement to numerical 
data and does not require an additional fitting parameter, as will be shown in Sec. |U 



4 Applications 

We study the direct regular-to-chaotic tunneling process in the systems introduced in 
Sec. |2] starting with the simplest example 2\o with a harmonic oscillator-like regular 
island. As further systems, we consider the map with a deformed regular region, the 
map T> rs with a regular stripe, and finally the standard map which shows a generic mixed 
phase space. 

The map 2\o has a particularly simple phase space structure with a harmonic oscillator- 
like regular island with elliptic invariant tori and constant rotation number, see the insets 
in Fig. |U Numerically, we determine tunneling rates by using absorbing boundary con- 
ditions at \q\ > 1/2. Analytically, for the approach derived in Sec. |3] Eq. fl27]) . we use 
the Hamiltonian of a harmonic oscillator as the fictitious integrable system U rEg . It is 
squeezed and tilted according to the linearized dynamics in the vicinity of the stable fixed 
point located at the center of the regular island. 

Figure |4] shows the prediction of Eq. ( 127]) compared to numerical data. We find 
excellent agreement over more than 10 orders of magnitude in 7. In the regime of large 
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1(T 

10 30 50 i/h eS 70 

Figure 7: Numerical data (dots) for m < 8 for the map with a harmonic oscillator-like 
island. Comparison with Eq. f l4"6]) (solid lines). The prediction of Refs. [El 05], Eq. fj49|) . 
for m = with a fitted prefactor is shown (dotted line). 



tunneling rates 7 small deviations occur which can be attributed to the influence of the 
chaotic sea on the regular states: These states are located on quantizing tori close to the 
border of the regular island and are affected by the regular-to-chaotic transition region. 
However, the deviations in this regime are smaller than a factor of two. 

In Fig. [7J we compare the results of the semiclassical prediction, Eq. (T4"6"j) . to the 
numerical data. Due to the approximations performed in the derivation of this formula 
stronger deviations are visible in the regime of large tunneling rates while the agreement 
in the semiclassical regime is still excellent. 

In Refs. [T5l H5] a prediction was derived for the tunneling rate of the regular ground 
state, 

r(a,4«) , Af ^ 

1 (a, 0) 

where T is the incomplete gamma function, a = A/h c s, and c is a constant. Eq. (|49l) can 
be approximated semiclassically (see Ref. [23]), leading to 

70 oc e h <=ff v . (50) 

Figure [7J shows the comparison of Eq. fj49l) (dotted line) to the numerical data for the 
map "Dho- Especially in the semiclassical limit deviations are visible. The factor 2 which 
appears in the exponent of Eq. is more accurate than the factor 3 — In 4 in Eq. flBU]) . 

In typical systems the rotation number of regular tori changes from the center of the 
regular region to its border which typically has a non-elliptic shape. Such a situation 
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Figure 8: Dynamical tunneling rates from a regular island to the chaotic sea for the map 
T>&. Numerical results (dots) and prediction following from Eq. ( 127)) (lines) vs l/h e g for 
quantum numbers m < 5. The insets show Husimi representations of the regular states 
m = and m = 5 at l//i e fl = 50. The prediction of Refs. [151 HH]> Eq. ( 14"9)) . for m = 
with a fitted prefactor is shown (dotted line). 



can be achieved using the map Here nonlinear resonances are small such that their 
influence on the tunneling process is expected only at large 1 / h e g, see the inset in Fig. [S] 
for its phase space. 

We determine the fictitious integrable system if reg by means of the Lie-transformation 
method described in Sec. 13.21 It is then quantized and its eigenfunctions are determined 
numerically. Fig. |H] shows a comparison of numerically determined tunneling rates (dots) 
to the prediction of Eq. (1271) (solid lines) yielding excellent agreement for tunneling rates 
7 > 10~ n . For smaller values of 7 deviations occur due to resonance-assisted tunneling 
which is caused by a small 10:1 resonance chain. Similar to the case of the harmonic 
oscillator-like island the fictitious integrable system U reg can be approximated by a kicked 
system U reg ~ UyU^ using V(q) = —rq 2 /2 + Rq 3 /3. Hence, Eqs. (135)) and (1361) can be 
evaluated giving similarly good agreement (not shown). The prediction of Eq. (149)) [T5"ll4"5] 
(dotted line) shows large deviations to the numerical data. 

We now consider the kicked system T) TS , for which the regular region consists of a stripe 
in phase space, see the inset in Fig. In Ref. [32] this map is used to study the evolution 
of a wave packet initially started in the regular region by means of complex paths. Also for 
this system one can predict tunneling rates by means of Eq. (1271) . The fictitious integrable 
system U TOg is determined by continuing the dynamics within \q\ < qd/(8ir) to the whole 
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Figure 9: Dynamical tunneling rates from a regular stripe to the chaotic sea for the map 
T> TS : We compare numerical results (dots) and the prediction following from Eq. ( 12 7p 
(lines) vs l/h e g for the quantum numbers |m| < 4. The inset shows the phase space 
of the system. The numerical data is obtained using an absorbing region at \q\ > 1/4 
(gray-shaded area of the inset). 



phase space. It is given as a kicked system, Eq. ([2]), defined by the functions 

^» - 4(»4-t)- < 5i > 

T'{p) = ~sm(2irp). (52) 

For sufficiently large smoothing parameter b, see Eq. ffTOj) . the dynamics of the mixed 
map inside the regular region is equivalent to that of the regular map. Thus Eq. 027|) 
can be applied and we compare its results to numerically determined data. Absorbing 
boundary conditions at \q\ > 1/2 lead to strong fluctuations of the numerically determined 
tunneling rates as a function of h e g, presumably due to dynamical localization. Choosing 
\q\ > 1/4 f° r the opening, which is closer to the regular stripe, we find smoothly decaying 
tunneling rates, see Fig. O The comparison with the theoretical prediction shows quite 
good agreement. Note, that due to the symmetry of the map there are always two 
regular states with comparable tunneling rates except for the ground state m = 0. These 
two states are located symmetrically around the center of the regular stripe. While the 
theoretical prediction, Eq. fl27|) . is identical for both of these states, the numerical results 
differ slightly due to the different chaotic dynamics in the vicinity of the left and right 
border of the regular region. 

The paradigmatic model of an area preserving map is the standard map, see Sec. 12.11 
For k between 2.5 and 3.0 one has a large generic regular island with a relatively small 
hierarchical region surrounded by a 4:1 resonance chain, see the inset in Fig. [TQl Absorbing 
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Figure 10: Tunneling rates for the standard map (k = 2.9) for m < 2. Prediction 
of Eq. (j27p (lines) and numerical results (dots), obtained using an absorbing region at 
\q\ > 1/4 (gray-shaded area of the inset). 



boundary conditions at \q\ > 1/2 lead to strong fluctuations of the numerically determined 
tunneling rates as a function of h c g, presumably caused by partial barriers. Choosing 
absorbing boundary conditions at \q\ > 1/4, which is closer to the island, we find smoothly 
decaying tunneling rates (dots in Fig. [TU]). Evaluating Eq. (|2T|) for k = 2.9 gives good 
agreement with these numerical data, see Fig. [TU] (solid lines). Here we determine H Teg 
using a method based on the frequency map analysis, as the Lie transformation is not able 
to reproduce the dynamics within the regular island of U, see Sec. 13 .21 With increasing 
order of the expansion series of H reg the tunneling rates following from Eq. fTSTl) diverge. 
Hence, for the predictions in Fig. [TU] we choose terms up to second order only. Note, 
that at such small order the accuracy of H rcg within the regular region of U is inferior 
compared to the examples discussed before. Hence, in Eq. (j2Tjl the state U\ij)™^) has small 
contributions of all purely regular states l^g) 111 ^ ne regular island. These contributions 
are removed by the application of the projector P c h. However, this projector depends on 
the number of regular states N reg , which grows with l/h cS . If N TCg increases by one, P TCg 
projets onto a larger region in phase space. This explains the steps of the theoretical 
prediction, Eq. (|2"7|) . visible in Fig. [TUJ 

5 Summary 

Dynamical tunneling plays an important role in many areas of physics. Therefore a de- 
tailed understanding and quantitative description is of great interest. In this text we 
have given an overview on determining direct regular-to-chaotic tunneling rates using the 
fictitious integrable system approach. The direct regular-to-chaotic tunneling mechanism 
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is valid in the regime where Planck's constant is large compared to additional structures 
in the regular island, such as nonlinear resonances. To include the effect of such reso- 
nances, resonance-assisted tunneling [2"0~] |2"2"] has to be considered in addition with the 
direct regular-to-chaotic tunneling contribution. Recently, the two mechanisms have been 
studied in a combined theory [23] , giving a full quantitative description of tunneling from 
a regular island, including resonances, into the chaotic region. 

The approach to determine direct regular-to- chaotic tunneling rates can also be gen- 
eralized to the case of billiard systems. It can be applied if the fictitious integrable 
system H reg is known, as e.g. for the annular or the mushroom billiard. In addition to 
H reg a description for the chaotic states of the mixed system H is needed, for which we 
employ random wave models [5] which account for the relevant boundary conditions of 
the billiard. For the mushroom billiard the fictitious integrable system is easily found 
as the semi-circle billiard and with this an analytical expression for the tunneling rates 
was derived in Ref. [18J. This result has been compared to experimental data obtained 
from microwave spectra of a mushroom billiard with adjustable stem height and good 
agreement was found. It was also shown that tunneling rates manifest themselves in ex- 
ponentially diverging localization lengths of nanowires with one-sided surface roughness 
in a perpendicular magnetic field [3S1 WI\ . 

The fictitious integrable system approach has been extended to open optical micro- 
cavities in Ref. [19]. In particular the annular microcavity was studied, which allows for 
unidirectional emission of light and shows modes of high quality factors simultaneously. 
This is desirable for most applications. In contrast to closed billiards the leakiness of 
the cavity has to be considered, which leads to the contribution Q^ T to the quality factor 
caused by the direct coupling of the regular mode to the continuum. Additionally, the con- 
tribution caused by dynamical tunneling Qdyn is relevant, where 1/Q — 1/Qdir + ^/Qd y n- 
Qdyn is directly related to the dynamical tunneling rates given by the theory using a fic- 
titious integrable system. The prediction for the quality factors Q has been compared to 
numerical data. Excellent agreement is found if no further phase-space structures exist 
in the chaotic sea. If additional structures appear, the numerical data show oscillations, 
which cannot be explained by the present theory. 

Future challenges include a completely semiclassical prediction of direct regular-to- 
chaotic tunneling rates in generic systems, the understanding of implications of additional 
phase space structures, and the extension to higher dimensional systems. 
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